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We compute exactly the low-energy single-electron Green's function, the impurity and electron 
self-energies, and the resistivity for the two-channel Anderson impurity model. These results are 
obtained by exploiting the boundary conformal field theory identified from the Bethe Ansatz so- 
lution of the model. Using that solution we can make contact with the parameters of the original 
Hamiltonian and provide the detailed crossover between the two integer valence limits. Our results 
generalize those obtained previously in the context of the two-channel Kondo model. 



I. INTRODUCTION 

Multi-channel impurity physics is a well-established 
route to non-Fermi-liquid behavior. Its appeal is there- 
fore understandable as an essential ingredient in many 
of the various theories that try to explain the unusual 
characteristics of numerous systems ranging from heavy 
fermionai to mesoscopic point contacts^ In this con- 
text, the two-channel Anderson impurity model occupies 
a central place. It was first introduced by Cox in an 
attempt to model the physics of certain U-based heavy 
fermions, 3 among which the compound UBei3 is a prime 
example. Following that line of thought, a dilute concen- 
tration of uranium in a ThBei3 matrix will correspond 
to a metallic system with a feeble concentration of two- 
channel impurity centers. Impurity corrections to the dif- 
ferent transport properties of that system should display 
fractional power laws indicative of non-Fermi-liquid be- 
havior. For instance, corrections to the resistivity would 
display \[T dependence at low temperatures^ and should 
constitute a particular experimentally accessible observ- 
able. 

Indeed, transport measurements are good candidates 
for experiments, since they can also be performed in 
mesoscopic systems for which bulk thermodynamic mea- 
surements, possible for heavy fermions and other materi- 
als, are ineffectual. In this context, the two-channel An- 
derson model was already used successfully as the start- 
ing point of Non-Crossing Approximation (NCA) calcu- 
lations to model the temperature dependence of the dif- 
ferential conductance of Cu point contactsA^ 

The two-channel Anderson impurity modeli&S de- 
scribes the interaction of three-dimensional electrons 
with a local impurity carrying both spin and quadrupolar 
degrees of freedom. These degrees of freedom correspond 
to the lowest-energy configurations of a uranium impurity 
in charge states U 4+ (5/ 2 ) and U 3+ (5/ 3 ). Taking into 



account spin-orbit coupling and crystal-field splitting in 
a cubic background one ends up considering the minimal 
scenario of a (quadrupolar) T3 doublet representation of 
the cubic group as the lowest-energy multiplet in U 4+ and 
a Kramers doublet Tg as the lowest-energy configuration 
in U 3+ . These levels hybridize with those electrons from 
the conduction band effectuating the transition from one 
doublet to the other. Starting from the full electron field 
\l/(x), one carries out an expansion in harmonics corre- 
sponding to cubic symmetry, retaining only ^& acr (r) which 
has the appropriate symmetry to couple to the impurity, 
with a = ± denoting the quadrupolar degrees of freedom 
and a =|, J, denoting magnetic (spin) degrees of freedom. 

We then proceed to write the field in terms of ID right 
(left) moving fields ip Raa (ipLaa-)' representing the incom- 
ing (outgoing) radial components of the 3D electron fields 
that couple to the impurity, 



2\/27rr 



-ikpr 



(1) 



with a "free electron" boundary condition ip Laa {0) 

ipR alT (Q) imposed at the origin r = 0. 10 
The Hamiltonian is then given by 



H = ifi 



bulk 1 -"ion 



hybr : 



(2) 



where 
#bulk 



dr 



^La( r )( i9 r)^W( r ) : - 
" : fL CT M(^r)Vw( r ) 



#ion = £sP a h + e q blb a 

ffhybr = V[^ Laa (0)by a + /> a Vw(0)] 



(3) 

(4) 
(5) 



The conduction electrons here hybridize with the impu- 
rity via a matrix element V jii with the impurity modeled 
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by a quadrupolar (magnetic) doublet of energy e q (e s ), 
created by a boson (fermion) operator b\ (/J) ^2 Strong 
Coulomb repulsion implies single occupancy of the local- 
ized levels: + b [ & b & = 1. The free part of the Hamil- 
tonian, -ffbulki defines a linearized spectrum around the 
Fermi level. The Fermi velocity is set to unity, with the 
resulting ID density of states p = l/(2n). Normal order- 
ing is taken with respect to the filled Fermi sea. 

The model in Eq. (0 was recently solved by two of 
us using a Bethe Ansatz*& A complete description of 
spectrum and thermodynamics was given, and it was 
found that at low temperatures the theory is attracted 
to a line of fixed points parametrized by the impurity 
charge valence n c (where n c — (flfa) measures the av- 
erage charge localized at the impurity site). Near in- 
tegral charge valence n c ~ 1 (n c ~ 0) a magnetic 
(quadrupolar) moment forms at intermediate tempera- 
tures. This moment is then screened by the conduction 
electrons as the temperature is lowered, leading to a zero- 
temperature entropy S™ lp = fcs In v2 and impurity spe- 
cific heat C" np ~ TlnT, typical of two-channel Kondo 
physics. In the mixed-valence regime one finds the same 
low-temperature behavior, but without the formation of 
a magnetic or quadrupolar moment at intermediate tem- 
peratures. 

In previous work we constructed the Boundary Con- 
formal Field Theory (BCFT) which describes the ap- 
proach to criticality J» The leading scaling operators were 
identified — including the exactly marginal operator that 
generates the line of fixed points — and all physical scales 
and BCFT parameters were determined explicitly via a 
numerical fit to the exact solution in Ref. 0]. This al- 
lowed us to go beyond the Bethe Ansatz approach and 
derive the critical exponents of the Fermi edge singu- 
larities caused by time-dependent hybridization between 
conduction electrons and impurity. Our results chal- 
lenged those obtained by more conventional, approximate 
schemes (see, e.g., Ref. 0). 

In the present work we take the BCFT formulation 
one step further and extract the exact space- and time- 
dependent single-electron Green's function of the model. 
This allows us to calculate the self-energies of the con- 
duction electrons and of the impurity, as well as the 
zero-temperature resistivity and leading temperature- 
dependent term. 

The analysis is most easily performed by generalizing 
that of the multi-channel Kondo model in Refs. [lfj and 
[r?|. In fact, the very structure of the Green's function, 
as well as that of the leading terms of the resistivity, can 
be read off directly from the corresponding result for the 
two-channel Kondo modeliiSiii The only essential new 
element in the analysis is how to properly introduce the 
scales and amplitudes that determine the influence from 
the magnetic and quadrupolar degrees of freedom as one 
moves away from the integer valence limits. 

In the n ext sec tion we combine results from Ref. fl4| 
and Rcfs. to obtain the single-electron Green's 

function of the model. In Sec. Ill we use this result to 



derive the self-energy of the impurity, and in Sec. IV we 
report on the calculation of the resistivity. Section V 
contains a summary and a discussion of our results. 



II. GREEN'S FUNCTION 

At sufficiently low energies (or large distances) a quan- 
tum impurity interacting with conformal (i.e., linear dis- 
persion) electrons can be represented by a conformally 
invariant boundary condition — this is the key idea of the 
BCFT formulation of a quantum impurity problem^ For 
the two-channel Anderson model, this boundary condi- 
tion [which supersedes the trivial "free-electron" bound- 
ary condition ip Laa (0) = ip Raa (0)], is most easily de- 
scribed via a "gluing condition" on the charge, spin, and 
flavor conformal towers that make up its spectrum (for 
details, see Ref. 14]). The dependence on the new 
boundary condition (alias the impurity) is picked up by 
the time-ordered left-right single-electron Green's func- 
tions 

G L R(T;r 1 ,r 2 ) = G RL (r;r 1 ,r 2 ) 

= {i> La «(T,ri)i>^ a (0,r 2 )) (6) 

(with t imaginary time) via the one-particle S'-matrix 

S {1) (u; F ) = e 2lS -C {1) (u JF ). (7) 

Here C^(ujf) is the amplitude for a single electron to 
scatter elastically off the impurity at the Fermi level ujf, 
and Sf — 8(lof) is the corresponding single-electron scat- 
tering phase shift. At large (mean) distances from the 
boundary, | r\ — r 2 | 3> a, with a some characteristic mi- 
croscopic scale, one finds tha1*i& 

5(i)(wf) 



GLi?(r;ri,r 2 ) 



(8) 



r + i(ri + r 2 ) 

Thus, at the level of the left-right Green's function, the 
boundary condition that emulates the presence of the 
impurity is coded by S^(ujf)- In contrast, the large- 
distance left-left (LL) and right-right (RR) Green's func- 
tions G mm (T;n,r 2 ) = (-^ma^ri^L^O,^)) with 
m = L, R are insensitive to the particular boundary con- 
dition imposed^ 



G LL (T;ri,r 2 ) = G* RR {r;ri,r 2 ) 
1 



r + i(ri - r 2 ) 



|ri-r 2 |>a. (9) 



Turning now to the the 3D electron field Green's func- 
tion G(T,r 1 ,r 2 ) = (# Q(T (t, ri)*J lCT (0, r 2 )) and expressing 
it in terms of the ID propagators in Eqs. © and © one 
obtains 

g(T,n,r 2 ) = e-^^- r ^G LL (T, n ,r 2 ) (10) 
+ e ik ^- r ^G RR (r, ri ,r 2 ) 
+ e- lk F^+^G LR (T, ril r 2 ) 
+ e ik ^ +r ^G RL (r,n,r 2 ). 
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Continuing the ID electron fields analytically to the 
full line — oo < r < oo [with ^u aa -(T,r) = 

i>L acr (T, —r)], averaging over impurity locations (thus 
restoring translational invariance), and exploiting a T- 
matrix formulation^ the Fourier-transformed Green's 
function in Eq. (|10|l can be cast in the standard form 



<7(w„, k) = 



1 



(11) 



with the self-energy S(w n ) given by 



E(w„) 



Ci(w F ) 



2irg F 



sgn(w„). (12) 



Here rii is the density of a dilute random distribution of 
impurities, qf is the 3D "free-electron" density of states 
at the Fermi level, and sgn(w„) is the sign function, with 
u) n = 2n(n + l/2)k B T, n = 0, ±1, ±2, . . ., Matsubara fre- 
quencies. In order to tie the self-energy £(u; n ) to the 
two-channel Anderson model in Eq. (|2J) we need to de- 
termine the 5-matrix in Eq. (|7|) that is associated with 
the hybridization interaction in Eq. (JSJ). The amplitude 
Ci\-\{u)f) can be determined in exact analogy with the 
two-channel Kondo problem in Ref. Within the 

BCFT formalism C\ (ojf) gets expressed as a certain com- 
bination of so-called "modular S matrices' — 



s: 



1 

71 



sm[n(2j f + l)(2j' f + l)/A] 



(13) 



with structure and allowed values of the quantum num- 
bers j f ,j' f = 0,1/2,1 determined by the SU(2) 2 Kac- 
Moody symmetry of the flavor sector— Specifically, 



Cx{w F ) 



pl/2o0 

CO q1/2' 
a l/2'-'0 



and it follows from Eq. (|13(l that 
d(w F )=0. 



(14) 



(15) 



Thus, the outgoing scattering state has no single-electron 
component after interaction with the impurity. This ex- 
treme non-Fermi-liquid behavior is the same as for the 
two-channel spin Kondo model^I (n c — 1 limit of the 
two-channel Anderson model) and is not modified as one 
moves into the mixed valence regime with n c ^ 0,1. As 
seen from Eq. (J7|), the impurity valence n c , connected to 
the phase shift dp via the Friedel-Langreth sum ruleS^ 



(16) 



could only influence the scattering if there were a finite 
single-electron cross section at the Fermi level. However, 
as C\{ujf) = independent of n c , this does not happen. 
[Note that the phase shift Sp is that of an electron with 
spin a and flavor index a, hence the unconventional fac- 
tor of 1/4 in Eq. lfTo|) . 14 ] 



To summarize the analysis thus far: The zero- 
temperature single-electron Green's function is given by 
Eq. (Uni , with 



1 



271-g j 



-sgn(w n ), 



(17) 



where ni and qf are defined after Eq. (|12fl . We should 
here stress that impurity-impurity interactions have been 
neglected in Eq. l(T7|) . As discussed in Ref. 0] this type 
of analysis is applicable only at temperatures high enough 
so that any remnant effects from interimpurity interac- 
tions are washed out (but low enough so that the the- 
ory is critical and the BCFT formulation remains valid) . 
Also note that the expressions for the ID propagators in 
Eqs. JSJ and @ pick up corrections when | r\ — r 2 \ < 
implying that the result in Eqs. lfTT|) and (|T7jl gets mod- 
ified when probing boundary correlations with large mo- 
menta. 

To obtain the leading finite-temperature and frequency 
corrections to Eq. i|17|) we need to consider the theory 
slightly off the line of boundary fixed points. The scaling 
Hamiltonian ^scaling that governs the critical behavior 
close to the fixed line was identified in Ref. 0] as 

HscaUng = H* + A c J(0) + X s O (s \0) + XfO^{0) 

+ sublcading terms . (18) 

Here H* is the critical Hamiltonian that represents Kbuik 
in Eq. (J2J, subject to the boundary condition that emu- 
lates the impurity terms H lon and ifhybr m Eqs. and 
(0, respectively. These terms, which break particle-hole 
symmetry, also give rise to the exactly marginal term 
A c J(0) in Eq. GBJ, with J(0) = E v : f aa ^ aa {Q) : 
being the charge current at the impurity site and with 
A c its conjugate scaling field. This is the operator that 
generates the line of stable fixed points. Off the fixed 
line the terms -f/ion an d ^hybr allow for additional ir- 
relevant boundary operators to enter the stage. Of 
these, A S 0( S >(O) and X f O^(0), both of scaling dimen- 
sion A = 3/2, are the leading ones. The spin boundary 
operator 0( S )(O) is the same operator that drives the crit- 
ical behavior in the two-channel spin Kondo problem and 
is obtained by contracting the spin-1 field 4)^(0) with 
the vector of SU(2) 2 raising operators : O^(0) = 
jLI • <t> (s) (0). The flavor boundary operator X f O^(0) 
has the same structure. In obvious notation: O^'(0) = 
J U l-<t> {f} (0). 

In the case of the two-channel (spin) Kondo problem 
the flavor operator is effectively suppressecUS: Of the two 
available energy scales, the bandwidth D and the Kondo 
temperature Tk (where Tk sets the scale for the crossover 
from weak coupling (high-temperature phase) to strong 
renormalizcd coupling (low-temperature phase)), only D 
enters the expression for the flavor scaling field A f . This 
is so since the Kondo temperature is dynamically gen- 
erated in the spin sector (as indicated by the infrared 
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divergences in perturbation theory) and hence cannot in- 
fluence the scaling of the flavor degrees of freedom. On 
dimensional grounds one concludes that A/ ~ 0(l/\/~D), 
whereas A s ~ 0(1/s/Tk). For a small Kondo coupling, 
call it A, Tk ~ £)exp(— 1/A) <C D, and the critical be- 
havior is therefore driven by (0) alone. As we showed 
in Ref. ^4[, the picture for the two-channel Anderson 
model is more involved. There are here two dynamically 
generated temperature scales T s j(e), both parametrized 
by e = e q — e s and thus varying with the position on 
the fixed line via the dependence of the impurity va- 
lence n c on The scaling fields X s j are parametrized 
accordingly^: 



-W - 



B 



(19) 



with B s j dimensionless constants. The pre cise depen- 
dence of A s j on e was extracted in Ref. [14j from a fit 
to the numerical solution of the thermodynamic Bethe 
Ansatz (TBA) equations of the model^ In the magnetic 
moment regime where e <C H — T [two-channel (spin) 
Kondo limit] T s <C Tf and A s dominates. As e increases 
T s and Tf approach each other and become equal when 
e = [i (maximal mixed valence with no moment forma- 
tion). Continuing along the fixed line, the two scales 
trade places, and eventually, at the quadrupolar critical 
end point (e » [i — Y) one finds that T s ^> Tf. It follows 
that in the two-channel Anderson model both boundary 
operators O^(0) and 0^(0) come into play, with their 
relative importance changing continuously as one moves 
along the fixed line. 

Going back to the scaling Hamiltonian in Eq. Q18JI. 
the effect of the exactly marginal term A c J(0) is eas- 
ily obtained via the observation that it samples the lo- 
cal charge at the impurity site, with the scaling field 
A c = — n c /4 measuring the impurity valence per spin 
and flavor degree of freedom»i4 By the Friedel-Langreth 
sum rule in Eq. (|l(j|l , the resulting shift of the charge con- 
tent of the critical bulk Hamiltonian H* , Q — > Q — n c , 
shows up as a phase shift Sp — 7m c /4 on the elec- 
trons that scatter off the impurity charge potential at 
r = 0. In other words, ipLaa — * exp(— iirn c /A)ipLaa 
and ipRaa —> exp(i7m c /4)?/>.R Q0 ., implying that the left- 
right propagators Glr = G* RL get phase shifted by 
25 p = irn c /2, as indicated in Eq. ||8J. 

To probe the effects from the spin and flavor boundary 
operators in Eq. (|18|) requires a perturbative approach. 
Passing to a Lagrangian formalism, the correction SS 
to the Euclidean fixed point action due to O^(0) and 
0( S )(O) in Eq. l(T%|l can be written as 



SS = X * 

k=f,8 



dT'j^l-(f)^(T',0), 



(20) 



with j3 — 1/ksT. To leading order in a perturbative 
expansion this leads to the following correction for the 



left-right propagator: 



S ai S„ tl 5G{T\r 1 ,r 2 ) = V" A fc / dr' 

k=f,s J ° 

x (V'w(^n)^ fe i ) -0 (&) «O)4 w (O,r 2 )) r . (21) 



The index T that appears in Eq. H21|l refers to the "finitc- 
T geometry" T + = {w = t + ir}, connected to the half- 
plane C + = {Imz > 0} used at zero temperature via the 
conformal mapping w — (/3/7r) arctan(z). 

The three-point functions in Eq. I|21[l are completely 
determined by conformal invariance up to multiplicative 
constants Nf and N s 



.16 



5G LR (T,ri,r 2 ) = i(XfN f + X s N s )e 



2M> 







isin|(r + i(n + r 2 ))Y 



X / dr' 

/o [sin 5(r' - r - in) sin |(r' + ir 2 )Y- 



(22) 



In the n c = (quadrupolar) limit where X s — + and 
Sp = 0^* the theory is invariant under charge conju- 
gation (particle- hole symmetry). Adapting an argument 
from Ref. 0] we may use this property to determine the 
phase of Nf, and together with an explicit calculation of 
| Nf | 2 one finds the value 



Nf = 3/V8. 



(23) 



The value of iV 2 is fixed via Eq. (|23J) by the vanishing of 
the four-point function^ 

(^L ff ( T i +* r i ) ( T i ~ ir i ) V>1 7M (T2 +ir 2 ) i> Rlli (r 2 -ir 2 )) 

]1/2 32 



^{Nj+N^inr,)^- 



n - r 2 



and one concludes that 



N s = -i 



V8' 



(25) 



with the negative sign in Eq. (|25(l following from the con- 
dition that the expression for SGlr in Eq. (|22|) collapses 
to that for the two-channel Kondo modeli& in the n c — > 1 
limit. The scaling fields A/ and A s are the same as those 
that parametrize the thermodynamics and can thus be 
fitted to the exact TBA solution of the model— (see the 
next section). With this fit SGlr in Eq. (|22|l will be com- 
pletely specified. Note that by time reversal invariance, 
SGrl = SG* LR . 

Turning to the chiral propagators in Eq. @ , it is easy 
to verify that the corrections to these from SS vanish 
identically: SGll (SGrr) is also given by the integral 
expression in Eq. I|22|) but with r 2 — ► — r 2 (n — > — n). 
All zeros of the denominator are located in the upper 
(lower) half plane, and the integration contour can be 
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deformed to Imr — > — oo (+00) without crossing any sin- 
gularity; hence 5G LL = SG RR = 0. 

The integral expression for SGlr = 5G* RL in Eq. lt2*2t 
differs from that for the two-channel (spin) Kondo model 
in Ref. only by having a prefactor (A/ — iX^)e™ na ^ 2 
instead of a single scaling field A s (ee A in Ref. [la]). (It 
follows trivially that the result for the chiral propaga- 
tors is the same for the two models.) From this point 
on we can therefore carry over the analysis intact from 
Ref. [3, at the end simply taking A (A/ -i\ s )e™ n ^ 2 . 
This gives for the finite-temperature and frequency self- 
energy^ 



rii sgn(w„) f 3 . x \ i—n f 2lT \ 2 

X /' I «» l/ 2 - W -V2 (1 _ „)V2 

Jo L 

- r- 2 (3/2) u- 1 ' 2 (1 - ur 3/2 j J , (26) 

where F (u) is the hypergeometric function 
F(3/2,3/2,l;«). 

To summarize this section: to leading order in temper- 
ature and frequency, the exact single-electron Matsubara 
Green's function Q (u> n ,k) of the two-channel Anderson 
model for a dilute distribution of impurities is given by 
Eq. JTUJ) with the self-energy S (uj n ) in Eq. These 
are bulk-electron quantities; in the next section we make 
an aside to discuss their connection with the impurity 
response. 
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FIG. 1: Imaginary part of the zero-temperature impurity self- 
energy as a function of frequency. Different curves correspond 
to different values of the impurity configurational energy split- 
ting (e). The frequencies are scaled with the Kondo temper- 
ature Tk- 



one following the same conventions as in Ref. |26j, the 
relation reads 



Z R (w) = n, 



V 2 
2wg F mp 



(30) 



Parametrizing the impurity Green's function with a spec- 
tral weight equal to 1/2, a hybridization amplitude 14 
T = npV 2 , and a self-energy E^ (w), we can extract 



III. IMPURITY SELF-ENERGY 

For the sake of simplicity, we will carry out this dis- 
cussion at zero temperature. By analytic continuation to 
real frequencies, 



lim E(w„) = E K (w), 



(27) 



one obtains from Eq. (|26[) an integral expression for the 
retarded electron self-energy Y, r (lo). By taking the T — > 
limit and approximating the integral as done in Ref. [l6l ] 
one finds tha'fcSi 



£t=o («>) 



with 



(28) 



12 



1 + — (A+a 2 ) (l-isgn(w)) m 5 

V7T 



A\ (A/, A s , n c ) = Xf cos (n c 7r/2) + X s sin (n c 7r/2) 
A2 (Xf,X s ,n c ) ee Xf sin (n c 7r/2) - A s cos (n c 7r/2) 

Writing down the relevant equations of motion, one can 
make a connection between the self-energy in Eq. I|28|) 
and the impurity Green's function^ Defining the latter 



iT 1 



[E*H] M . (31) 



The resulting formula for the impurity self-energy in- 
herits from E fl (u>) a range of validity for \u>\ -C Tk, 
where Tr- ee 4r/7r 2 e _ ^' £ '^ r is the two-channel Kondo 
temperature^ 

In Fig. ^ we show the imaginary part of the retarded 
impurity self-energy as a function of frequency (scaled 
with the Kondo temperature) for different values of the 
energy splitting e. Only the curves for positive val- 
ues of this parameter are shown in the figure, since 
for negative values one simply has to use the relation 
ImE^pfw,- e) = ImS^f-w, e) to obtain the corre- 
sponding curves. These results are in fair agreement with 
the ones obtained recently using the Numerical Rcnor- 
malization Group (NRG) mrthodiSi Additionally, notice 
that the figure required knowledge of A\ and A2 as func- 
tions of e; we shall comment on this in the next section. 



IV. RESISTIVITY 

Given the retarded electron self-energy T, r (oj), defined 
by Eqs. I|26|) and l|27|l . one can readily obtain the resis- 
tivity p{T) of the model to leading order in temperature. 
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Adapting the argument in Appendix C of Ref. 0], the 
assumption that the impurity-electron interaction is well 
described using the two-channel s-wave decomposition in 
Eq. 0} implies that the resistivity can be expressed di- 
rectly in terms of lmE R (uj). From Eqs. (|2"BT and itTTjl one 
obtains 



ImS R (w) = - 



1 



Al(A/, A s ,n c )x 



x / du 
lo 



2ng F { ' ' V/?„ 
cos(/3w(lnw)/27r)u^(l - u)?F(u) 

- T- 2 (3/2)u-i(l-uY 

+ A 2 (X fl X s ,n c )x 



<l du sin (/3oj (In it) /2tt) u *(1-u) 2 F{u) 
lo 



(32) 



We have here used the property that 

i/l^l/^sgn^) -> cos(f^lnu) - isin(^lnu) 

under the analytic continuation iuj n — > to + i0 + . 

As shown in Ref. |25j for this class of problems, vertex 
corrections to the resistivity involve s-wave correlations 
which vanish identically. It follows that the Kubo for- 
mula for the resistivity contains only the quasiparticle 
lifetime r, with no weighting over large-angle scattering 
processes. With two local orbital channels (a = ±) of 
charge conduction, the Kubo formula thus readsi^ 



P-\T) 



4e 2 
3m P 



d 3 k 



dn F (e k ) 
de k 



k 2 r(e k ). (33) 



Here e and m e are the electron charge and mass, re- 
spectively, iif(ft) is the Fermi distribution function, 
and T(efe) is the lifetime of a quasiparticle of energy 
£fc = k 2 /2m e , 

T{e k ) = -±[lmX R (e k )]-\ (34) 
Combining Eqs. JS^J, lO, and it follows that 
4:TT(eg F v F ) 2 „/7r\ 5 f°° dx 



du 



, 4cosh 2 (x/2) 
Ai{Xf, A s ,n c ) 

cos(a;(lnu)/27r)w~5(i _ F(u) 

- r- 2 (3/2)iH(i- u )- 

+ „4 2 (A/,A s ,n c ) 



; / sin(x(lnu)/27r)u ^(l — u) 2 F(u) 



(35) 



where x = e^/^T '. We have here used that [dn^/de^] in 
Eq. (|33|l rapidly goes to zero away from the Fermi level, 



0.1 



T7T K 



FIG. 2: Resistivity vs temperature curves for different values 
of the microscopic Hamiltonian parameter e. Positive and 
negative values of e fall on top of each other; from left to 
right the curves correspond to |e| ~ 0.0, 0.5, 0.8, 1.1, 1.4, 1.8, 
and 2.1. The curves corresponding to the last two values of 
e fall on top of each other in the scale of the figure, which 
illustrates the collapse into a universal curve as the system 
goes away from mixed valence. 



allowing us to approximate the momenta that appear in 
the integral by k F (= m e v F ). (Note that previously v F 
was set to unity. We still use units where h = 1.) Carry- 
ing out the integrals over x and subsequently over u (this 
second integral has to be done numerically, but to ma- 
chine accuracy the result is found to be a rational number 
and expected to be exact) one finally obtains, for the 
low-temperature resistivity, 



p(T)=p(0) 



with 



P(0) 



-) Ai(Xf,X s ,n c ) 



4ir {eg F v F y 



(36) 



(37) 



Since both A/ and X s on the one hancU^ and n c on the 
other hancU^ are known to be functions of e, the ex- 
pression for the leading low-temperature dependence is 
related to the original impurity Hamiltonian via this sin- 
gle parameter. This behavior is illustrated in Fig. [21 
While we have plotted the curves over the full interval 
< T < Tk, we should alert the reader that our results 
are exact only in the scaling regime T <C Tk- The ex- 
plicit connection between A/, s and e was obtained from a 
fit of the low-temperature thermodynamics to the results 
of the TBA solution^ But the thermodynamics involves 
only the squares of the scaling holds and the sign remains 
therefore undetermined. However, in the limit of integer 
valence we have 



lim Ai(X f (e),X s (e),n c (e)) 

e/r—*±oo 



A/,, (e) , (38) 



7 



and in that limit the model maps onto the weak-coupling 
two-channel Kondo model (Jk < J* , with J* the Kondo 
fixed point under renormalization to low energies) for 
which the sign of A/ iS (e) is known to be negative (it 
is expected to reverse sign for Jk > J*) Ml Since the 
combined BCFT and TBA analysis indicates that the 
scaling fields vary continuously and do not change sign, 
we conclude that the coefficient Ai is always negative. 
Hence, the resistivity is a monotonically decreasing func- 
tion of temperature for all values of e. In the case of n c , 
the Bethe Ansatz solution provides directly an expression 
that relates it to 

With these considerations the values of Ai^s) are 
completely determined, which allowed us in the previ- 
ous section to plot the impurity self-energy. It is inter- 
esting to point out that a positive sign for the scaling 
fields (A/ jS ) will hamper the comparison of that plot with 
the NRG results. Moreover, a positive sign will spoil 
the causal properties of the self-energy and is therefore 
unphysical for the Anderson model (even in the mixed- 
valence regime). 

Away from e sa 0, we can use Eq. (|19|l and our results 
of Ref. |14| to derive an approximate expression for the 
resistivity that highlights its scaling properties, 



p(T)/p(0) = l 



cos (n c 7r/2) + sin (n c 7r/2)] 




Here T K = min {T f ,T s } is the BCFT Kondo scale and 
the prefactors are in correspondence with the precise def- 
inition of this scale as given in our previous work. It is 
important to remark again that the asymptotic expan- 
sion for small temperatures is valid only for T <C Tk- 
For consistency, the same definition of the Kondo scaled 
was used for both plots in Figs. 1 and 2 (the prefactors 
in the definition are of course not universal and depend 
on the particular convention). 

It is here interesting to discuss the experimentally 
measured low-T resistivity of the thoriated UBei3 com- 
pound, mentioned in the Introduction, which shows a 
Vt behavior bu t wit h a positive coefficient. As dis- 
cussed in Refs. |4l27j . this would imply that in the 
(single-impurity) Kondo model framework this system 
exhibits a strong electron- impurity coupling (Jk ></*). 
It was speculated that such a regime was achievable near 
mixed valence in the context of the (single-impurity) two- 
channel Anderson model (see Ref. [7j for a review). Our 
results, however, do not support those ideas. Perhaps, 
since Ui-^Th^Beia with x = 0.1 is far from the di- 
lute limit, lattice effects might play a role in reversing 



the sign of that coefficient. While later measurements 
on the same compound have confirmed the vT behav- 
ior with a positive coefficientr& on the other hand, a 
y/T scaling of the resistivity with a negative coefficient 
has been recently observed in a different uranium-based 
heavy-fermion material, Sci^^U^Pda (but only for large 
dopings, x m 0.35j2£ when a single-impurity description 
of the low-temperature physics is not always expected). 



V. SUMMARY AND DISCUSSION 

In a previous article 1 ^ we presented the technical de- 
tails of the BCFT solution of the two-channel Ander- 
son model and discussed its asymptotic low-temperature 
thermodynamics. In that context we were able to make 
the connection with the full solution obtained using the 
Thermodynamic Bethe-Ansatz formalism and could ex- 
plicitly match the scaling fields with the microscopic pa- 
rameters of the lattice Hamiltonian. Here we completed 
the task by calculating dynamical- and transport-related 
quantities (the single-electron Green's function, the elec- 
tron and impurity self-energies, and the resistivity). Us- 
ing our previous results from Ref. ^4[, these quantities 
are parametrized directly in terms of the energy dif- 
ference between impurity configurations in the original 
Hamiltonian (e). We have shown, in particular, how our 
analytic expression for the impurity self-energy captures 
the low-frequency behavior in agreement with the results 
of other nonperturbative techniques like Wilson's numer- 
ical renormalization group method^ As we mentioned 
in the Introduction, having reliable access to transport 
quantities is of crucial importance for the comparison and 
interpretation of experiments that continue to seek indis- 
putable realizations of multi-channel Kondo physics. The 
work presented here furthers our understanding of two- 
channel Kondo physics in mixed- valent scenarios, thus 
widening the range of possible candidate systems for ex- 
perimental realizations. 
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